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1 Introduction 



The Hubbard model is an important model in condensed matter physics. Lieb and 
Wu [Q] diagonalized the 1-D Hubbard model in terms of the coordinate Bethe Ansatz. 
The existence of the Bethe Ansatz equations implies the integrability of the model. 
In order to prove it from the QISM, Shastry f|, §] proposed a coupled 6-vertex model 
and constructed the transfer matrix related to the 1-D Hubbard model through the 
R matrix and the L operator. Olmedilla, Wadati and Akutsu [§], ||, || obtained the 
super symmetric R matrix and L operator by applying the Jordan- Wigner transforma- 
tion to the L operator in || |J. In ||, the eigenvalue of the transfer matrix for the 1-D 
Hubbard model was conjectured. Using the coordinate Bethe Ansatz method, Bariev 
derived the eigenvalue of the diagonal-to-diagonal transfer matrix of the coupled 
6-vertex model. Recently, we []TDj] found the eigenvalues of the transfer matrices of the 
1-D Hubbard model and the coupled 6-vertex model with twisted boundary condition 
by using the Analytic Bethe Ansatz method. The eigenvalue of the transfer matrix of 
the 1-D Hubbard model was also studied by Ramos and Martins (ll[ from the viewpoint 
of the Algebraic Bethe Ansatz method. 

The R matrix of the 1-D Hubbard model has the quite different property. The 
-R(/ii,/i 2 ) matrix of the XXZ model depends on the difference of the two spectral 
parameters /ii and /x 2 , where /Zj is defined for the 2-th auxiliary space; If [i\ — /i 2 = 

— /i' 2 , then i?(/ix,/i 2 ) = R([i[, /4) • For the 1-D Hubbard model, however, this is not 
the case. Furthermore, the R matrix and the L operator of the 1-D Hubbard model 
are different; the R matrix depends on two spectral parameters, while the L depends 
on one. As printed out in f||], the proof of the Yang-Baxter equation RRR = RRR 
is independent from that of RLL = LLR. Shiroishi and Wadati proved the Yang- 
Baxter equation RRR = RRR and obtained a new hamiltonian with periodic boundary 
condition. In this paper, we discuss the diagonalization of the inhomogeneous 2-D 
model with twisted boundary. By applying the Analytic Bethe Ansatz approach, we 
find the eigenvalue of the transfer matrix and the Bethe Ansatz equations. Under the 
uniform limit (homogeneous case), the logarithmic derivative of the eigenvalue of the 
transfer matrix for the 2-D model gives the energy spectra of the derived hamiltonian 
with twisted boundary. If the twisted angles are zero, the hamiltonian is reduced to 
that given in @. 
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2 Model 

The R matrix related to the Hubbard model is given in [|| 

8 sinh(/i! + /i 2 ) 
c/ sm2(/ii + /i 2 ) cos(/xi - /x 2 ) 

x {cos(/ii + /x 2 ) cosh(/ii - h 2 )L1 2 (fi 1 - [i 2 )Ll 2 (ni - H2) ^ 

+ cos(/ii - /x 2 ) sinh(/ii - /i 2 )L^ 2 (/ii + /x 2 )cr 2 L[ 2 (/xi + ^2)^} , 

where <r and r are Pauli matrices. The /Xj are spectral parameters and hj describing 
the interaction strength. They are controlled by sinh2/tj = (U sin 4. The L a for 
a = a, r are defined by 



L\ 2 {n) = tu 4 (/x) + w 3 (h)<t^<t 
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L[ 2 (/x) = tw 4 (/i) +w 3 (h)t*t% + r 1 + r 2 + r x r 2 + , (2) 

where 

w 4 (/x) + u> 3 (/x) = cos(xx), 

w 4 (/x) - w 3 (fj) = sin(/x). (3) 
This R matrix satisfies the Yang-Baxter equation 

^2l(/"l, ^2)^31(^1, A*3)-R32(A*2, ^3) = -R31O1, ^3)^32(^2, ^3)^21 (/^l, A*2)- (4) 

The Yang-Baxter equation (4) implies the existence of a general 2-D inhomogeneous 
model with the transfer matrix 

{/Xj}) = tr g T(n, {/ij}) = tr g R Lg (n, /x L ) ■ ■ ■ Ri g (n, /ii). 

The logarithmic derivative of £(/x, {/Xj}) at /x = /ij = /xo gives a 1-D quantum system 
with periodic boundary condition 0. Here, we consider the twisted transfer matrix 

t(MM) = r 11 (/x,{/x i })e-^ + T 22 (/x,{/x i }) e -^ 

+r 3 3(/x, W)e^ + T 44 (/i, {/x,})^ (5) 

where and are free parameters (twisted angles). It is easy to prove that the 
twisted transfer matrix is the generating function of the infinite number of the conserved 
quantities. 

With the initial condition of R matrix 

Rmgi^O, [I'm = /^o) = Pmg = Pmg^mgi (6) 



being arbitrary parameter, the hamiltonian related to the twisted transfer matrix 
(5) is given by the logarithmic derivative 

L-1 L-1 

H = E ( a m+l a m + °m+l°m) + E ( T m+l T m + T m+l T m) 
m=l m=l 

+ Sin(2/i ) Om + l°m - <Wl^m) } ( 7 ) 

x {cos 2 {ho)t^ - sin 2 (/i )r^ +1 + sin(2/i )(T+ +1 T- - r m+1 r+)} 

+ 4cosh(2/> ) {^oK-sm^K 

+ sin(2/i )(e i ^ + ^a L "(Tf - e^^V^fff)} 
x {cos 2 (/i )r£ - sin 2 (/i )rf + sin(2/i )(e i( ^ ) r L F rf - e^-^rZ rf)} . 

This hamiltonian has four free parameters, the c7 and /xq denote the interaction strength. 
The and ?/> describe the twisted boundary. This hamiltonian will reduce into the 1-D 
coupled XY model under /j = 0. Generally, the hamiltonian (7) contains the interac- 
tion between the charge sector and spin sector. This is due to that the R matrix does 
not depends on the difference of two spectral parameters. Notice that this hamiltonian 
will recover the one in [^] at <fi — ip — 0. 



3 Diagonalization 



In ||10|| , we obtained the eigenvalues of the transfer matrices related to the 1-D Hubbard 
model and the coupled twisted XY model by using the Analytic Bethe Ansatz method. 
Here we want to find the eigenvalue of the transfer matrix (5) by taking use of the 
same idea. Let us define the total reference state to be the state with all spins down 
(spin a and spin r). The matrix T(/i, {//./}) on the reference state takes the form 



T(/i, {[ij})\vac >- 



A^n) 

* A 2 (/i) 

* A 2 (» 

* * * A 4 ([i) 



\vac >, 



(8) 
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where * stands for the no- vanishing terms and 

L 

3=1 
L 

L 

e h-hj cos (jj^ cos(p,j) — e hj ~ h sin(/x) sin(/ij 



s(/w) 



cos 2 (/i) — sin 2 (/i J ) 
-e h ~ hj cos(/i) cos(/i,) - e hil ~ h sin(/i) sin(//_ 7 -), 
pg(/x, /ij) = sin(/i — fij) cosh(/i — ^) — sin(/x + /x^) sinh(/i — hj), 
Pi(/x, /Xj) = cos(/i — /ij) cosh(/i — hj) + cos(/x + /Xj) sinh(/i — /ij). (9) 

Using the explicit expression of the R matrix, one can easily find the eigenvalues of the 
states with N r-spin (or cr-spin) flipping from the reference state. For the 1-D Hubbard 
model, N corresponds to the number of electrons. After a long but direct calculation, 
we arrive at 

+AMe-»Il , Pl \ M ( ) (10) 



where 



e h hj cos(/x) cos(uj) — e hj sin(/x) sin(z/ J ) 



cos z {j2) — sin [y 



3> 

p 4 (/x, Vj) = e h] ~ h cos(/x) cos(i^) + e h ' hj sin(/x) sin(i^), 
, . e~ 2h] cos(/x) sin(/x) — e~ 2h cos(vj) cos(uj) 

P&KPiVj) = — — r , 

cos 2 (/x) — sm \Vj) 
p 3 (/x, z/j) = e ft ~ Aj ' sin(/x) cos(^) — e hj ~ h cos(/x) sin(p,). 

where the notation smh(2hj) = U sin(2i/j)/4 has been used. 
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In accordance with the hypothesis of the Analytic Bethe Ansatz approach [jTU 
now seeks for a more general form 

atv(^) = A 4 (^n^4 



one 
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=1 fj) 



m=l 



where g^O^, A m ) and ^(/i, A m ) are some undetermined functions. Here iV is the total 
number of cr up-spins and r up-spins, M the number of r up-spins. Because the R 
matrix has no crossing symmetry, these functions can not be fixed by using the standard 
Analytic Bethe Ansatz method. Let us consider some properties of these functions. 
The eigenvalue A(/i) as the analytic function of [i should have vanishing residues at 
the poles fi = uj and // = i/j, which sets up the relation 

M M 

II 9s( u j^m) = II 92(pj,Xm) (13) 
m=l m=l 

where e _2hj cot(Pj) = e 2h cot(^). Second, from the special case of A at N = 2, M — 1, 
we know that the functions #2 and #3 are the rational and have some poles. In terms 
of the news variables k defined by e lk = —e~ 2h cot(/i), g<i and g% can be written as 

92 ^ A = • ■ m — \~TTTu ' 93 A = • ■ 77s — r— 7777 14 
1 sm(/cj — A + U/4 1 sm(/cj — A + (7/4 

where P2 and P3 are analytic function of \i. The equation (13), together with the 
asymptotic behavior of A(/x), fixes completely P2 and P3. Thus, we obtain the final 
result 

AM = AMe^fi^^- 

, j / u# TT P^ v i) tt isin(fc) -A m -?7/4 
2l/iJ M^(^^)iii zsin(A;)-A m + ^/4 
+A 2 (//)e-^ TT PioO^j) tt »sin(fc)-A m + 3Z7/4 1 j 
2 J=i PiO", fj) - Ps{p, Vj) iii i sm(k) - A m + C//4 

+^(Ai)e-^n , ^fr^i v 



:i2) 
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the Bethe Ansatz equations are 

" ismjk^-Xn + U/A t^M+lj*/, tt A TO - A m - U/2 

^[ismik^-Xn-U/A { ] r ll i A„-A m + C//2' 



(16) 



where e tkj = —e~ 2hj cot(uj). Equations (15) and (16) are the eigenvalue and the Bethe 
Ansatz equations for the 2-D inhomogeneous model. If the solutions of equation (16) 
are given, we can find the eigenvalues of A. Furthermore, we can calculate the ther- 
modynamics of the system. 

Now, let us return to consider the energy spectra of the hamiltonian (7). Taking 
the logarithmic derivative of A(/i) at /i — /i and letting /ij — /i , we get 

E = ^ 2 ^ UL ^-2Ecos(.,) + ^E 

4 2cosh(2/i ) K 3> 2fy 

| cosh~ 1 (2/i )[sinh(2/i — 2hj) — cos(2z/j) cosh(2/i — 2hj)] 



. cot(2/xo) sin(2z/j) + sinh(2/i — 2hj) — cos(2uj) cosh(2/i — 2hj) 
[sinh(2/i — 2hj) — cos(2i/,-) cosh(2/i — 2hj)\ 
[cot(2/x ) sin(2i/,-) + sinh(2/i — 2^-) — cos(2z/j) cosh(2/i — 2hj)] 
[cosh(2/ij) + cos(2z/j) smh(2hj)] 



(17) 



x 



sin(2i/,-) 

sinh(2/ij) + cos(2z/j) cosh(2^) 



cot(2/i ) sin(2i/j) + sinh(2/i — 2/ij) — cos(2z/,) cosh(2/i — 2/^) J 

4 conclusion 

We have found the eigenvalue and the Bethe Ansatz equations for the 2-D inhomoge- 
neous twisted model by making use of the Analytic Bethe Ansatz approach. We have 
also obtained the energy spectra (17) for the related 1-D quantum hamiltonian (7). 
In equation (17), the first three terms are similar to that of the 1-D Hubbard model. 
The terms in the curved bracket reflect the contribution due to the interaction between 
the spin and charge sectors. In the R((/,i, 1^2) matrix, the spectral shift does not keep 
the invariance of R, i.e. R(/ii + 5, ^2 + S) 7^ R(/ii,{i 2 )- Thus the interaction between 
the two sectors can not removed by using the spectral shift. The uniform shift of z//s 
can not give the uniform shift of k/s due to the nonlinear relation between i/^'s and 
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k/s. This is very different from the case in which the R(fii, H2) just depends on the 
difference of the two spectral parameters. 

In fact, the eigenvalues of the 1-D Hubbard model and the coupled XY model can 
be considered as the special case of equation (17). The L operator of the 1-D Hubbard 
model is obtained from the Ripx, ^2) by setting /i 2 = in 0, [|. This means the 
uniform parameter /j, being zero. At this case, the equations (16) and (17) recover the 
results given in [|i"0|j . 

Based upon the equations (15)-(17), one can calculate the finite-size correction and 
get the conformal scales. Notice that the twisted angles <fi and ip can be interpreted as 
the external vector potentials coupled to the system, one can consider the conductivity 
of the system as done in 1-D Hubbard model and 1-D t- J model |12, [13| . 
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